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Self-consistent solutions to nonlinear spinor field equations 
in General Relativity are studied for the case of Bianchi type- 
I space-time. It has been shown that introduction of A-term 
in the Lagrangian generates oscillations of the Bianchi type-I 
model. 
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1. INTRODUCTION 

The quantum field theory in curved space-time has 
been a matter of great interest in recent years because 
of its applications to cosmology and astrophysics. The 
evidence of existence of strong gravitational fields in our 
Universe led to the study of the quantum effects of ma- 
terial fields in external classical gravitational field. After 
the appearance of Parker's paper on scalar fields Q and 
spin-i fields, Q several authors have studied this subject. 
Although the Universe seems homogenous and isotropic 
at present, there are no observational data guarantying 
the isotropy in the era prior to the recombination. In fact, 
there are theoretical arguments that sustain the existence 
of an anisotropic phase that approaches an isotropic one 
j|. Interest in studying Klein-Gordon and Dirac equa- 
tions in anisotropic models has increased since Hu and 
Parker Q have shown that the creation of scalar particles 
in anisotropic backgrounds can dissipate the anisotropy 
as the Universe expands. 

A Bianchi type-I (B-I) Universe, being the straightfor- 
ward generalization of the flat Robertson- Walker (RW) 
Universe, is one of the simplest models of an anisotropic 
Universe that describes a homogenous and spatially flat 
Universe. Unlike the RW Universe which has the same 
scale factor for each of the three spatial directions, a 
B-I Universe has a different scale factor in each direc- 
tion, thereby introducing an anisotropy to the system. It 
moreover has the agreeable property that near the singu- 



larity it behaves like a Kasner Universe, even in the pres- 
ence of matter, and consequently falls within the general 
analysis of the singularity given by Belinskii et al ||] . Also 
in a Universe filled with matter for p — je, 7 < 1, it 
has been shown that any initial anisotropy in a B-I Uni- 
verse quickly dies away and a B-I Universe eventually 
evolves into a RW Universe Since the present-day 
Universe is surprisingly isotropic, this feature of the B- 
I Universe makes it a prime candidate for studying the 
possible effects of an anisotropy in the early Universe 
on present-day observations. In light of the importance 
of mentioned above, several authors have studied linear 
spinor field equations |plp| and the behavior of gravita- 
tional waves (GW's) [pppd| in a B-I Universe. Nonlinear 
spinor field (NLSF) in external cosmological gravitational 
field was first studied by G.N. Shikin in 1991 @. This 
study was extended by us for the more general case where 
we consider the nonlinear term as an arbitrary function 
of all possible invariants generated from spinor bilinear 
forms. In that paper we also studied the possibility of 
elimination of initial singularity especially for the Kas- 
ner Universe . For few years we studied the behavior 
of self-consistent NLSF in a B-I Universe [HJl5| both m 
presence of perfect fluid and without it that was followed 
by the Refs., Jl|Hl| where we studied the self-consistent 
system of interacting spinor and scalar fields. The pur- 
pose of the paper is to study the role of the cosmological 
constant (A) in the Lagrangian which together with New- 
ton's gravitational constant (G) is considered as the fun- 
damental constants in Einstein's theory of gravity E^] . 



2. FUNDAMENTAL EQUATIONS AND 
GENERAL SOLUTIONS 

The Lagrangian for the self-consistent system of spinor 
and gravitational fields can be written as 



L = 



R + 2A 
2k 



rmptp - 



(2.1) 
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with R- being the scalar curvature, k- being the Ein- 
stein's gravitational constant. The nonlinear term 
describes the self-interaction of a spinor field and in 
this particular case is chosen as an arbitrary function 
of 5 = $^, i.e. L N = F{S). 

Variation of (2.1) with respect to spinor field gives non- 
linear Dirac equations 



«7 M V^ - $V = 0, 
iV M ^ + = 0, 



(2.2a) 
(2.2b) 



with $ = m — dF/dS, whereas variation of (2.1) with re- 
spect to metric tensor gives the Einstein's field equa- 
tion 



■MS 



(2.3) 



where B% is the Ricci tensor; R = g iil 'R ilv is the Ricci 
scalar; and Tjf is the energy-momentum tensor of matter 
field defined as 



rpp 

t 1 



4 
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(2.4) 



j sp with respect to (2.2) takes the form 



L, 



OF OF 

dip dip 



F. 



(2.5) 



In fl2.2| ) and (2.4) V M denotes the covariant derivative 
of spinor, having the form [^0 21 



dip 



(2.6) 



where r^(x) are spinor affine connection matrices. 7 ma- 
trices in the above equations are connected with the flat 
space-time Dirac matrices 7 in the following way 

9^{x) = el(x)el(x)r) ab , ^(x) = e a ^(x)j a , 

where r\ a b — diag(l, —1, —1, —1) and e° is a set of tetrad 
4- vectors. 

Let us now choose Bianchi type-I space-time in the 
form I2I 



ds 2 = dt 2 - a 2 dx 2 - b 2 dy 2 - c 2 dz 2 



(2.7) 



with a, b and c being the functions of t only. 

For the space-time (2.7) Einstein equations ( |2.3| ) now 
read 
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-A, 


(2.8a) 


%-kGTI 


-A, 


(2.8b) 


8irGTi 


-A, 


(2.8c) 


8vrGT 


-A, 


(2.8d) 



where point means differentiation with respect to t. 
From 



(2.9) 



one finds 



To = 0, Tt 



AQrr, 



(2.10) 



r 2 = ^W7 2 7 ,r 3 = ^)7 3 7 . 



Flat space-time matrices 7 we will choose in the form, 
given in |2^| . 

We will study the space-independent solutions to 
spinor field equation (7) so that ip = ip{t). Setting 



abc ■ 



-9 



in this case equation (2.2a) we write as 



'V'( §- t + ^ ]"-<!><' = ()■ 



(2.11) 



(2.12) 



Furth er putting Vj(t) = y/ripj(t), j = 1,2,3,4, from 
(|2.12 ) one deduces the following system of equations: 



Vi + i<f>Vi = 0, 
Va ~ i$V 3 = 0, 

V 2 + i$V 2 = 0, 
Vi - i^V A = 0. 



From ( 2.13 ) one easily obtains 



Wi,2 = —pre 



1 ^3,4 in 
V3,4 = —;=e , 



(2.13a) 
(2.13b) 

(2.14) 



where f2 = [ Q dt and Cj are the constants of integration. 

From ( 2.12Q we will also find the equation for bilinear 
spinor form S = ipip: 



with the solution 



S+ -S = 0, 

T 



(2.15) 



(2.16) 



where Co is the constant of integration. As one can see 
the constants Co and Cj are connected with each other 
as 

r< _ r<2 1 ri2 n2 r<2 



Putting ( 2.16 ) into (2.4), we obtain the following ex- 
pressions for the components of the energy-momentum 
tensor 
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rpO 
1 



mS-F(S), Tl=T, 



dF 



(2.17) 



From (2.17) in view of (2.16) it is obvious that for the 
linear spinor field 

T ° = mS = mCo/r, = T| = T 3 3 = 0. (2.18) 



The sign of Co is defined from the requirement of posi- 
tivity of energy density Tq of linear spinor field. Hence, 
from ( 2.18 1 ) em erges Co > 0. In view of (2.17) from (2.S) 
we derive ||13H15| 



o(t) = (#^3)1/3^/3 exp [ (2Xl + x 3 ) 0], (2.19a) 
b(t) = ( J D 1 - 1 D 3 ) 1 / 3 r 1 / 3 exp[-(Xi - X 3 ) <&], (2.19b) 
c(i) = (L>i^)- 1/3 T 1/3 exp[-(Xi + 2X 3 ) (2.19c) 

where we define § = J [l/3r(t)]dt. 

As one sees, both spinor field and metric functions are 
in some functional dependent on r. Let us defi ne th e 
equation for t. Summation of Einstein equations (2.< 
(|2.8bj),([2~8c]) and Q2.8dD multiplied by 3 gives 



3kt 



The right-hand-side of 
namely 

3k/ _ dF„ 



T °j - 3Ar (2.20) 
2C|) is a function of r only, 



2Ft) -3At:= T{t) 



the solution to this equation is well-known for any arbi- 
trary function T(t) urn and can be written in quadrature 



dr 



= t + t 



2jT(r)dT 



(2.21) 



where to is some constant that can be set zero. Given the 
explicit form of the nonlinear term F(S) from (2.21) one 
finds the concrete solution for t. Thus the initial systems 
of Einstein and Dirac equations have been completely 
integrated. 



3. ANALYSIS OF THE RESULTS 

In this section we shall analyze the general results ob- 
tained in the previous section for concrete nonlinear term. 

Let us consider the nonlinear term as a power function 
of r, precisely F(S) — XS n with A being the coupling 
constant and n > 1. Inserting F(S) into (2.20) one ob- 
tains 



3kCq 



m + A(n-2)-| zr 



- 3A 7 



(3.1) 



The first integral of (|3.]|) has the form 



3kCq 



n — l /n — 2 



TO 



3Ar 2 



(3.2) 



with To being some positive constant. Finally, we obtain 

r («-2)/2 dr 

= = ±V3t. 

^kCqTot 11 - 2 + KC mT n - 1 - At" - ACJ 

(3.3) 

Depending on the sign of A and A we have the following 
pictures. 

case 1. A = — e 2 < 0, A > 0. In this case for n > 2 
and t — > 00 we find 



T(t) 



,V3et 



(3.4) 



Thus we see that the asymptotic behavior of r doe s not 
depend on n and defined by A - term. From (2.19) it is 
obvious t hat the asymptotic isotropization takes place. 

From ( |3.3| ) it also follows that t cannot be zero at any 
moment, since the intigrant turns out to be imaginary 
as t approaches to zero. Thus the solution obtained is 
a nonsingular one thanks to the nonlinear term in the 
Dirac equation and asymptotically isotropic. 

Let us go back to the energy density of spinor field. 
From 



rpO _ 

1 o — 



mC AC,? 



follows that at 



T n-i < XC o 



(3.5) 



(3.6) 



the energy density of spinor field becomes negative, which 
means that the absence of initial singularity in the consid- 
ered cosmological solution appears to be consistent with 
the violation of the dominant energy condition in the 
Hawking-Penrose theorem p2[ , since in this case 



(3.7) 



T} = Tl = T 3 = A( " ~ > 0. 



Consider the linear case with A = 0. Then from ( |3.3| ) 
follows 



Tit) = -e 
V ' 4 



As one sees 



, -V3et ( K 2 m 2 Cl _ kCoTo \ 

+ v 4e 4 e 2 ; 



KmCo 



lim r 

t — >oo 



1 

I 6 



V3et 



that coincides with (3.4). from ( p.q ) follows 
1/ KmCo\ 2 kC t 



(3.8) 



(3.9) 



(3.10) 



that means, r(0) is defined by the relation between the 
constants. 
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case 2. A > and A > 0. For n > 2 ( ^3| ) admits 
only nonsingular oscillating solutions, since r > and 
bound from above. Consider the case with n = 4 and for 



simplicity set to = 0. Then from (3.2) one gets 



r(t) 



1 

Ia 



kC t 



J o 'o 



1/2 



(3.11) 



For a spinor field with A > and A > and n = 10 
a perspective view of r is shown in FIG. 1. Period for 
massive field (to ^ 0) is greater than that for massless 
one (to = 0). The initial value of r has been taken to 
be unit, i.e., r(0) = 1. For r(0) = 10~ 2 for example, r 
increases tremendously at initial steps (from 10~ 2 to 10 13 
in first step) and amplitude in that case is 10 16 , whereas 
the value of n (order of nonlinearity) defines the period 
(the more is n the less is period). 




FIG. 1. Perspective view of r showing the initially nonsin- 
gular and oscillating behavior of the solutions 

case 3. A < and A < 0. The solution is singular at 
initial moment, that is 



[^-Z\n 2 C%/At] 2,n (3.12) 
and at t — > oo asymptotic isotropization takes place since 



lim t 

t->0 



lim r w e 

t — >oo 



/3At 



(3.13) 



case 4. A > a nd A < 0. Solution is initially singular 
and coincides with ( 3.12 ) and bound from the above, i.e., 
oscillating, since 



(3.14) 



lim r w sinv 3At. 

t— »oo 



4. CONCLUSION 

Within the framework of the simplest nonlinear model 
of spinor field it has been shown that the A term plays 
very important role in Bianchi-I cosmology. In particu- 
lar, it invokes oscillations in the model which is not the 
case when A term remain absent. Growing interest in 
studying the role A term by present day physicists of 
various discipline witnesses its fundamental value. For 
details on time depending A term one may consult pq ] 
and references therein. 
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